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Abstract. Modified versions of two “standard” pulsar search techniques 
are presented that allow large-scale searches for pulsations in long dura¬ 
tion high-energy data sets using relatively modest amounts of computer 
time. For small numbers of photons (N p ^ lot 10 4 ), optimized brute-force 
epoch folding searches are preferred. For larger numbers of photons, ad¬ 
vanced Fourier domain acceleration searches are used. Using these tech¬ 
niques, my collaborators and I have searched Chandra observations of the 
Cas A supernova remnant (SNR) point source and the isolated neutron 
star RX J1856.5—3754 for pulsations, and confirmed the 65.6 ms pulsar 
in the 3C 58 SNR during a blind search of archival RXTE data. 


1. Introduction 


Sensitive searches for pulsars in long duration high-energy (i.e. X-ray and 7 - 
ray) observations are difficult. For relatively short observations ( T 0 b s 20 ks) 
all known isolated pulsars can be assumed to have a constant spin period during 
the observation. In this limiting case, searches for relatively bright pulsars are 
usually implemented using simple epoch folding or a standard Fourier analysis 
(e.g. Leahy et al., 1983| ; van der Klis, 1989|) . When T 0 i, s is long, however, the 
spin-down of young high-energy pulsars becomes important. The corresponding 
change in frequency during an observation can make these pulsars invisible to 
traditional search techniques. 

Another difficulty involved in searching very long observations (i.e. many 
days) is the extraordinarily large number of independent trials required for 
a blind search. For epoch folding searches that cover a wide range of both 
frequency / and frequency derivative /, the computational complexity is oc 
NphotTob s , or oc T° bs for a source of constant flux. Advanced Fourier analy¬ 
ses of very long time series (i.e. (2 — 20 ) x 10 s points where the photons have 


been binned) have been carried out by a few groups (e.g. Middleditch et al. 
2000|; |Chandler et al., 2001) and are, in general, more efficient due to the use 


of the FFT (computational complexity oc T° bs log T a b s ). The difficulties involved 
in computing many FFTs of this size, however, have relegated such analyses to 
those with access to large-scale super-computing. 

The techniques that I describe here allow large-scale blind searches of very 
long high-energy observations to be carried out on modest (albeit modern) work¬ 
stations or inexpensive workstation clusters. 
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2. Optimized Epoch Folding (N p h 0 t 10 4 ) 


Epoch folding involves placing each event (i.e. photon) into one of a relatively 
small number of pulse phase bins (usually 2—20) by calculating the exact pulse 
phase of the event based on a trial pulsar ephemeris (e.g. / and /) and an 
arbitrary starting phase of the profile (i.e. phase zero). For each trial profile, 
one of many possible (but preferably computationally “cheap”) statistical tests 
is performed on the profile to determine if it shows a pulsed signal or not. 

For blind searches, a large range of / and / must be covered at high reso¬ 
lution in both the / and / directions. In order to maximize sensitivity to low 
duty-cycle (i.e. narrow) pulse profiles, one must oversample the Independent 
Fourier Spacing (1 /T obs ) in frequency and the equivalent independent spacing 
in / (~ 6 /T° bs ) by a factor of 4 or more. 

More advanced searches also attempt to match the complexity of the ex¬ 
pected pulse shape with an appropriate number of pulse phase bins at the op¬ 
timal starting phase of the profile. For example, sinusoidal profiles only require 
2 properly phased bins for good sensitivity, but narrow pulse profiles are best 
detected if a single phase bin has the same width as the pulse and is centered 
on it. Unfortunately, searching over various numbers of phase bins and starting 
phases greatly increases the computer time required in an already costly search. 

I have made modifications to the Bayesian method of epoch folding pre¬ 
sented by Gregory fe Loredo (1992 , 1996| ) that incorporates everything described 
above in an efficient manner. These modifications allow brute force searches of 
“reasonable” high-energy data sets in a few weeks or less on a fast workstation 
(i.e. N phot ~ 1000 photons, T obs £ 30 days, 0 < / & 500 Hz, |/| ^ 10 -9 Hzs -1 ). 

The calculations required to search over various numbers of phase bins and 
starting phases have been reduced dramatically from the exact method described 
in Gregory &; Loredo (1992). The basic idea is to initially fold the data into 60 
or 120 phase bins and then combine various numbers of these bins (i.e. 2, 3, 4, 
5, 6, 10, 12, 15, 20 bins) in an efficient manner at each of the 60 or 120 starting 
phases in order to maximize the signal-to-noise. 

For very long duration observations, huge savings in CPU time can be had 
by searching the photons in a hierarchical manner and accepting some loss in 
sensitivity. Since for a search of a continuously sampled observation of a source 
with constant flux the required number of operations is oc T^ b , searching only 
a short but continuous part of the data takes much less time than searching the 
full observation (the number of photons and required resolution in / and / are 
reduced significantly). Candidates above a very low threshold (say the top few 
percent) in this initial search are then examined using twice as many photons (or 
approximately twice the duration) at a correspondingly finer resolution. This 
process continues until all the photons are examined and the full resolution of 
the data is reached around the candidate. 

My collaborators and I have recently applied these techniques to brute- 
force searches of Chandra observations of two very interesting X-ray sources: the 
point source found near the center of the Cas A SNR ( Tananbaum, 1999|) an d 
the isolated neutron s tar RX J1856.5—3754 ( Walter, Wolk, Sz Neuhauser, 1996 ). 
Murray et al. (2002a|) searched a 50 ks HRC-S observation (OBSID 1857) from 
the Cas A point source and reported no convincing candidates. More recently, 
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we have searched the follow-up 50 ks observation (OBSID 1038) of Cas A and 
could neither find any convincing new candidates nor confirm any of the low 
significance candidates from the earlier search (manuscript in preparation). A 
similarly unsuccessful search for pulsations (and a corresponding 99% confidence 
level upper limit to pulsations of 4.5%) from the 450 ks Director’s Discretionary 
Time observation of RX J1856.5—3754 was reported by |Ransom, Gaensler, & 
Slanc (20021) . 


3. Fourier-Domain Acceleration Searches ( N p hot ^ 10 3 ) 


When the observed number of photons is too large to allow efficient epoch folding 
searches to be conducted, the photons can be binned into a time series and 
searched using acceleration searches (i.e. / is constant during the observation) 
instead. Time domain acceleration searches have been used by numerous groups 
for some time (e.g. Middleditch & Kristian, 1984) and involve either re-sampling 
or re-binning the original uncorrected data to account for a non-zero /. Each 
time series corresponding to a different / trial then requires a long FFT and a 
subsequent search for pulsati ons. Unfortunately, F FTs of very long time series 
are difficult to compute (e.g. Chandler ct ah, 2001 ). 

By using Fourier domain acceleration search techniques (Ransom, 2001), 
only a single FFT of the original time series is required. All / trials are computed 
by correlating complex template responses with the raw Fourier amplitudes from 
the original FFT (i.e. matched filtering in the Fourier domain). The single long 
FFT that is required can be easily computed using out-of-core memory FFT 
techniques on even modest workstations. The subsequent Fourier domain accel¬ 
eration search is very efficient since the matched filtering is inherently memory 
local and easily computed using short FFTs. 

My collaborators and I have used a Fourier domain acceleration search 
code that includes Fourier interpolation (to minimize the effects of “scalloping”, 
e.g. van der Klis, 19891 ) and harmonic summing (to improve sensitivity to low 
duty-cycle pulsations) to successfully search numerous long radio observations 
of globular clusters and X-ray observations of potential X-ray pulsars. New 
pulsars detected using these methods include the b inary millisecond radi o pulsar 
PSR J1807—2459 in globular cluster NGC65 44 (Ransom et al., 2001| ) and a 


probably isolated 4.714 ms pulsar in Terzan 5 (Ransom, 2001). 

Recently, [Murray et al. (2002b| ) used these techniques to confirm the new 

65.6 ms X-ray pulsar discovered with Chandra in the 3C 58 SNR and measure 
its spin-down characteristics. We conducted a blind Fourier domain acceleration 
search on an archival 36 ks RXTE observation (OBSID 20259-02-01-00 with 

20.6 ks on-source) and uncovered a ~ 6.4 a candidate as a sum of 8 harmonics 
(the strongest having only ~ 13 times the local mean power — showing the 
need for a harmonic sum). The pulse profile from the RXTE observation is 
shown in Figure 1. By comparing the measured pulse periods from the RXTE 
and Chandra data we were able to measure the period derivative (1.93 x 10~ 13 ), 
spin-down E (~ 2.7 x 10 37 ergs s -1 ), surface magnetic field strength (~ 3.6 x 
10 12 gauss), and characteristic age (~ 5380 y), showing that it is a surprisingly 
“normal” young pulsar. 
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2 Pulses of Best Profile3C58_RXTE_20259_2-20kev.toa Search Information 



Candidate: 65.66ms_Cand 

Telescope: RXTE 

Epoch, = N/A 

Epoch = 50721.78649603110 

T sample = 0.000657 

Data Folded = 34344192 

Data StdDev = 0.01166 

Profile Avg = 4671 
Profile StdDev = 68.34 



P, (ms) = N/A P. (ms) = 65-65923001(16) 

P'top* (s/s)L= N/A P* (s/s) = 3.39(57)x10-' 5 

P ".opp (s/s 2 ) = N/A P" bor , (s/s 2 ) = 0.0(1.6)x10 -17 

P orb (s> y - N/A e = N/A 

a,sin(i)/c (s) = N/A (rad) = N/A 

T peri = N/A 



Freq - 15.230149 (Hz) 

-4x1 CT* -2x10~ 5 0 2x1 CT* 4x10 -5 



” Educed 

- 

f: 


Period - 65.65923765 (ms) 


Figure 1 The 2—20kev RXTE profile of the young 65.6 ms X-ray pulsar J0205+6449 
recently discovered in supernova remnant 3C 58 (see §3). 
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